Compactons and semi-compactons in the extreme baby 

Skyrme model 

J.M. Speight* 
School of Mathematics, University of Leeds 
Leeds LS2 9JT, England 



Abstract 

The static baby Skyrme model is investigated in the extreme hmit where the energy 
functional contains only the potential and Skyrme terms, but not the Dirichlet energy 
term. It is shown that the model with potential V = |(1 + fs)^ possesses solutions with 
extremely unusual localization properties, which we call semi-compactons. These min- 
imize energy in the degree 1 homotopy class, have support contained in a semi-infinite 
rectangular strip, and decay along the length of the strip as By gluing together 

several semi-compactons, it is shown that every homotopy class has linearly stable solu- 
tions of arbitrarily high, but quantized, energy. For various other choices of potential, 
compactons are constructed with support in a closed disk, or in a closed annulus. In 
the latter case, one can construct higher winding compactons, and complicated super- 
positions in which several closed string-like compactons are nested within one another. 
The constructions make heavy use of the invariance of the model under area-preserving 
diffeomorphisms, and of a topological lower energy bound, both of which are established 
in a general geometric setting. All the solutions presented are classical, that is, they are 
(at least) twice continuously differentiable and satisfy the Euler-Lagrange equation of 
the model everywhere. 



1 Introduction 

Solitons are stable, spatially localized solutions of nonlinear field theories. Ordinarily, "spa- 
tially localized" means that the field ip{t, x) approaches some constant vacuum value (fo asymp- 
totically as \x\ oo, usually exponentially in |a;| (e.g. KdV and sine-Gordon solitons, abelian 
Higgs vortices) sometimes as a power law (e.g. sigma model lumps, instantons). How- 
ever, there are some systems where the solitons' spatial localization is much more severe: 
<f{x) = i^Q exactly outside some compact region of space. Since they have compact support, 
these solitons are called compactons in the literature. They were first discovered in general- 
ized KdV equations [8] , then in nonlinear Klein-Gordon models with W-shaped potentials [3] , 
that is, potentials with two degenerate vacua at each of which the potential has a V-shaped 
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singularity. All these compactons live in one (spatial) dimension. Moving to two dimensions, 
compactons have been constructed for the baby-Skyrme model, with energy density 

S = ^A9j(/7 • di^ + X dj^) ■ {di(p x djip) + V{(p) (1.1) 

again in the case where V{(p) has a V-shaped singularity at the vacuum, for example V = 
-^T^^^ [1]. Here ip : M."^ ^ S"^ G M.^ and A is a positive constant. Since V is singular, it is not 
surprising that the compactons are singular, and this makes their interpretation somewhat 
problematic. In what sense, precisely, are they solutions of the model? 

In both [3] and [1], the singularity of V can be interpreted as V having infinite second 
derivative at the vacuum, so that the "mesons" of the theory (propagating small perturbations 
about the vacuum) have infinite mass. An alternative mechanism to give the mesons infinite 
mass, without making V singular, is to take the hmit A — )■ in the above model [4, 2]. This 
limit, variously called the pure, restricted or (as in this paper) extreme baby-Skyrme model, 
has the interesting property of being invariant under all area-preserving diffeomorphisms of the 
spatial plane, and is claimed to have applications in condensed matter physics [4]. For various 
choices of (continuously differentiable) potential V, it has been found to support compactons. 
One problem with all the compactons found in [4] , and some of those found in [2] , is that they 
are not even once continuously differentiable, so again it is not clear in what sense they are 
solutions of the model. Certainly they are not classical solutions of the field equation, which 
is a second order nonlinear PDE. They may be solutions in the weaker sense that they locally 
extremize the energy functional, but to make precise sense of this is rather technical: given 
that the fields themselves are not C^, what should be the allowed space of variations (usually 
taken to be with compact support)? Since the compactons in [2] saturate a topological 
lower energy bound, it is likely that a precise formulation of their status as solutions is possible. 
The compactons in [4] are more problematic. 

In this paper, we begin by analyzing the extreme baby Skyrme model in a rather general 
geometric setting, taking physical space to be any orientable two-manifold M and target 
space to be any compact Riemann surface A^" (so the primary case of interest is M = 
and N = S"^ C M."^). The potential will be taken to be of the form V = \U'^ where C/ is a 
non-negative function on A^ with isolated zeros. By a solution of the model we will strictly 
mean a twice continuously differentiable map M ^ N satisfying the Euler-Langrange equation 
for E everywhere. In this setting, we prove a topological lower energy bound, saturated by 
solutions of a first order "Bogomol'nyi" equation. Solutions of this equation have a natural 
interpretation as area-preserving maps from (part of) M to (almost all) A^, with respect to 
a deformed area form on A^ (determined by U). We show that all Bogomol'nyi solutions are 
solutions of the field equation and conversely (on M = M^) that all solutions of the field 
equation are (piecewise) Bogomol'nyi. The bound is a generalization of various special cases 
discovered previously [2, 5, 7, 9], and our main contribution here is to place these results within 
a geometric framework, and give a geometric interpretation of the Bogomol'nyi equation. 

We then consider the specific case M = M^, A^ = S*^, U = l+^ps in detail. By exploiting the 
model's symmetry under area-preserving diffeomorphisms, we construct a degree 1 solution of 
this model with extremely unusual localization properties, which we call a semi-compacton. 
This solution is constant outside a semi-infinite rectangular strip, decays hke x~^°^^ along the 
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length of the strip, and minimizes energy within its homotopy class. By gluing several (anti-) 
semi-compactons together, we show that in every homotopy class the model has solutions 
of arbitrarily high (but quantized) energy, all of which are at least marginally stable (in 
particular, they are not saddle points). Wc also prove that the critical set of any solution 
of the model can have no bounded connected components so, in particular, solutions can 
never have isolated critical points. We compare our results with those of Adam et al [2], who 
construct exponentially localized fields in this model, clarifying precisely when their fields are 
solutions in the strong sense used here. 

We go on to consider various cases where U is not but V — ^U'^ still is, which is 
enough for the critical parts of the general theory to survive, giving a necessary condition 
on U for the existence of compactons. In the case U = (1 + (^3)", | < a < 1, we give a 
geometric construction of the compactons obtained in [2] , and show how their key quahtative 
features (e.g. energy and area) can be found without solving any equations. In the case 
U — (1 + (^3)^(1 — 933)'', where | < a,^ < 1, we construct annular (or closed string-hke) 
compacton solutions which minimize energy in their homotopy class, generalizing results in 
[2] (which correspond to the degenerate case where the annulus is a punctured disk). 

In the final section, we consider the model on a compact domain, with U — 0, showing 
that generically the model has no nontrivial solutions at all. We conclude by suggesting some 
interesting open questions concerning the dynamics of semi-compactons. 



2 The Bogomol'nyi argument 

It is convenient to place the extreme baby Skyrme model within a more general geometric 
framework. The model has a single scalar field (p : M ^ N where (M, g) is an oriented 
Riemannian two-manifold, representing physical space, and {N, h, J) is a compact Riemann 
surface (with metric h and almost complex structure J), the target space. Denote by uj = 
h{J-,-) the Kahler form (or area form) on A^. Let U : N ^ [0, 00) be a function with 
isolated zeros, the vacua of the theory. In this section, we will take M either to be compact, 
or to be Euclidean M^. In the latter case (which is of most direct physical interest) we impose 
the boundary condition 

(p{x) (po ^ C/~-^(0) as |x| ^ 00 (2.1) 

sufficiently fast that Jj^ip*u! converges. Throughout, cp is assumed to be at least C^. The 
energy functional of the model is 

E{^) - If \^*uj\^yo\M + l [ {UocpfvolM (2.2) 

= ly*ujf + l\\Uoipf = E4{cf) + Eoicp) (2.3) 

where voIm denotes the volume form on M, and we have introduced the notation || • || for the 
norm of a function, or form, on M. It will be convenient to denote the associated inner 
product by (•, •), so for forms a,p E J1*'(M), 

{a,/3)^ / aA*/3, ||q;||^ = (a, a), (2.4) 
Jm 
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where * : Q^^M) ^(M) is the Hodge map. To obtain the usual baby Skyrme model, 

one chooses N — S"^ C R^, the unit sphere, with the induced metric and with almost complex 
structure J : T^S^ T^S^, X ^cpxX,so that uj{X, Y) ^ {cp x X) ■ Y ^ cp ■ {X x Y). In 
this case, with respect to any oriented local coordinate system {xi,X2) on M, 

We begin by establishing a topological lower energy bound for E = E4 + Eq of Bogomol'nyi 
type. The argument has been discovered in particular cases by several authors [2, 5, 7, 9], and 
our aim here is to place these results in a general geometric framework. 

Proposition 1 For all cp : M ^ N, 

E{cp) > ±{U) [ 

JM 

where (U) is the average value of U on N, with equality if and only if 
Proof: Clearly 

< h\*^*uj^Uo^f = E^ I ^*{Uuj). (2.6) 
2 Jm 

By dimensions, Ucu is a. closed 2-form and, since H'^{N) = M, there exists a constant a G M 
and a e (N) such that 

Uuj^au + da. (2.7) 

Then 

{00, Uuj) = a\\ujf + {uj, da) = aYo\{N) + {5uj, a) = aYo\{N) (2.8) 
since u is coclosed. But 

(w, Uu) ^ I Uu^ {U)\o\{N) (2.9) 

J N 

where {U) denotes the average value oiU : N ^"K. The result immediately follows. □ 

We remark that, since ui is closed, ip*ijj is a homotopy invariant of </?. In the case of 
most interest, N — S'^, the bound becomes 

E{(p) >4:7r{U)\n\ (2.10) 

where n G Z is the degree of (/?. 

The Bogomol'nyi equation Lp*u = *U oip has an interesting geometric interpretation which 
we will use frequently in later sections. Let A^o = ^^~^(0) C iV, the set of vacua of the model, 
and N' — N\No, the target space with the vacua removed. We can equip N' with a deformed 
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area form fl = cu/U. Note tliat tliis area form blows up as one approaches Nq, the boundary 
of N'. Given a map ip : M ^ N, denote by its critical set, that is 

= {x e M : rank(d(^^) < 2}. (2.11) 

At any x e M<^, {(p*uj)x — 0, since we can always evaluate this 2-form on a basis of vectors 
one of which is in kerd</7j;. Hence, any solution of the Bogomol'nyi equation maps into Nq 
(sends critical points to vacua), and on M' — M\M^ satisfies 

= = *1 = voIm- (2.12) 

U o if 

That is (as observed for a special case in [7]): 

Reiricirk 2 Bogomol'nyi solutions are area preserving maps from (M', voIm) to {N' , Q). 

Note that, as usual, the Bogomol'nyi equation is a nonlinear first order PDE for (/?. This 
is in contrast to the Euler-Lagrangc equation for which is second order. In analogy with 
harmonic map theory, it is convenient to make the following definition. 

Definition 3 The tension field oi cp : M ^ N is 

r{(p) = -Jd(p^S{(p*uj) + (UgiadU) o (p. 

Here S — — *d* : Q*'(M) — >■ Qp~^(M), the coderivative adjoint to d, and ji denotes the metric 
isomorphism T*M — )■ TM induced by g. Note that T{ip) is a section of ip~^TN, the vector 
bundle over M with fibre T^(^x)N over x e M. We will also consistently denote the form 
*(p*uj by : M ^ M, so 

ip*io = F^voIm- 

Given a variation Lpt of Lp, with infinitesimal generator X = dt(pt\t=o G r{(p~^TN) a straight- 
forward calculation [10] shows that 

= {X,T{ip))= f h{X,T{^))vo\M. (2.13) 
Hence, the Euler-Lagrange equation is 

T{ip) = 0. (2.14) 

Any solution of the Bogomol'nyi equation 

F^^±Uoip (2.15) 

minimizes energy in its homotopy class, so must satisfy the field equation (2.14) by the fun- 
damental lemma of the calculus of variations. It is reassuring to verify this fact directly. The 
key observation is contained in the following lemma. 
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Lemma 4 Let (p : M ^ N and X be a vector field on M. Then 

h(d^X, t((^)) = -^d(F^ - {U o ^f)X. 

Proof: One sees that 

^S(p*uj = -^*d*(p*uj = dF^ = -Jm gradF^ (2.16) 

where Jm is the almost complex structure induced by the orientation on M. Hence 

h{dipX, r(9?)) = h{dipX, Jd^pJu grad F^p + [U grad U) o y?) 

= -^*u{X, Jm grad F^) + {UdU){d^X) 

= -F^g{X, grad F^) + {U o cp)g{X, grad(C/ o cp)) 

= -l9{X,gre.d{F'^-{Uocpy)) = -^d{F'^-{Uo^f)X. (2.17) 

□ 

We remark that this Lemma remains true under the weaker assumption that V = ^U^ is 
(rather than U itself). One replaces t/gradt/ and UdU by gradF and dV throughout the 
proof. 

Proposition 5 Let if : M ^ N satisfy (either of the) Bogomol'nyi equation(s), (p*u) = 
± * U o everywhere. Then ip satisfies the field equation T{(p) = 0. 

Proof: By assumption F^ — [U o ip)'^ is constant on M, so by Lemma 4 we have that 
h{dp^X,T{ip){x)) = for all x G M and aU X e T^M. It follows that T{ip){x) = at aU 
regular points of </?, since dip^iT^M) — T^N at such x. It remains to show that </? satisfies (2.14) 
on its critical set. So, let x be a critical point of (/? (meaning rankd(/?a; < 2). Then p)*uix — 
so F^p{x) = 0, and hence (p{x) G U~^{0). But [/ > 0, so (pi^x) is a minimum of U, and hence 
{gmdU){p{x)) = = df/^(^) also. Hence (gradF^)(a;) = tl(dF^)x = ±tldf/^(^)d(^j; = 0, and 
one sees from equation (2.16) that 5ip*u) = at a;. Hence (p satisfies (2.14) at x. □ 

So solutions of the Bogomol'nyi equation automatically satisfy the field equation, as usual. 
In a general field theory of Bogomol'nyi type, there is no reason why solutions of the field 
equation should necessarily satisfy the Bogomol'nyi equation. Remarkably, we will show that, 
on M = M^, all solutions of the field equation satisfy one or other of the the Bogomol'nyi 
equations at each point. 

Proposition 6 Let(p :M? ^ N satisfy the field equation (2.14) <J'i^d boundary condition (2.1). 
Then 

Fi = {u o 

everywhere 
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Proof: Since t((/9) = 0, we see from Lemma 4 that — [U o is constant on M?. But 
^{x) ifo E U^^{0) as \x\ — )■ oo by (2.1), so F^{x) and U{ip{x)) — )■ as \x\ — >■ oo. 
Hence, this constant is 0. □ 

Since its proof uses only Lemma 4, Proposition 6 extends immediately to the weaker 
case that V = ^U^ is C^. By contrast, the proof of Proposition 5 makes essential use of 
the differentiability of U, so does not extend to this weaker case. By the support of a map 
(p : M ^ N we mean the closure of ip~^{N'), that is 

supp(^ = closure({a; e M : (p{x) ^ Nq}). (2-18) 

It follows immediately from Proposition 6 that, for a solution (/? on M = R^, suppc/? = 

closure(M\M^). 

As we will see later, it is possible, for suitable U, to construct solutions to (2.14) by gluing 
together maps with ip*uj = *U o cp and (p*u! = — * U o (p in different regions of M. So it does 
not follow from Proposition 6 that all solutions of the theory are global energy minimizers. In 
particular, the vacuum sector can contain infinitely many static solutions, of arbitrarily high 
energy. Remarkably, we will see that these "lump-antilump" superpositions are actually local 
minima of E, not saddle points. A key property which we will exploit in the construction of 
these exotic multilumps is the invariance of the model under area preserving diffeomorphisms 
of {M,g). Once again, this property has been observed previously in specific cases by many 
authors [4, 7]. 

Proposition 7 Let ip : M ^ N and £/ : M' M be an area preserving diffeomorphism. 
Then E{(po£/) = E{(p). 

Proof: By assumption, =e/*volM = voIm'- Let ip = ip o . Then 

i\)*uj = s^*{p)*u) = ^*(F^vo1m) = {F^ o ^)vo1m'. (2.19) 

Hence 

EM - \ I (F^o ^)'yo\M' [ ^*{F^yo\M) ^ \ I F>o1m = E^ (2.20) 

^ JM' ^ JM' ^ JM 

since =2/ : M' — > M is a diffeomorphism. Similarly, 

^ JM' ^ JM' 

- \ l {U op)\o\M = Eo{^). (2.21) 

^ JM 

□ 

It follows immediately that r((/9 0£/) = t{p) os^ so po^/ satisfies the field equation (2.14) 
if and only if p> does. Since F^oj^ = F^o ^ we also verify immediately that p) o satisfies 
the Bogomol'nyi equation (2.15) if and only if p) does. 
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3 Semi-compactons 

In this section we restrict attention to the case M = R^, N — S"^ and 

U{ip) = l + ip3, (3.1) 

though the constructions below clearly generaUze to any U which is C^, non- negative and 
has a single non-degenerate zero. It is straightforward [2] to find a degree 1 solution of the 
Bogomol'nyi equation (2.15) within the hedgehog ansatz 

ip(r, e) = {s/l - z{rY cos 9, ^1 - z{tY sin z(r)) (3.2) 

where (r, 9) are polar coordinates on M^, z : [0, oo) — )■ M, z{0) = 1 and z{oo) = —1. In terms of 
cylindrical coordinates Z — (p^ and © = arg(9?i +iip2), the Kahler form on 5"^ is a; = d© AdZ, 
and the ansatz (3.2) is Q — 9, Z — z{r). Hence, the Bogomol'nyi equation (2.15) becomes 

--^1 + z (3.3) 
r 

whose solution, with the required boundary data, is 

^(r) = -l + 2e-'-'/'. (3.4) 

Note that this solution has faster than exponential decay, and is smooth everywhere, including 
at the origin. To check this, define the (globally) analytic function 

n=0 ^ ' 

and note that 

^/l^^^{^^2r{l-r\{r^)f^/^. (3.6) 
Since q{Q) ^ 0, \fq{^ is analytic in a neighbourhood of 0, so 

Vl--2(r)2 = rg(a:2 + y2) (3.7) 

where Q is analytic on a neighbourhood of 0. It follows that 

^{x,y) = (xQ(x2 + y^),t/Q(x2 + t/2),_i + 2e-(-'+^')/2) (3.8) 

is smooth at (0, 0). Clearly, (?7) = 1, so this unit lump solution has energy E — An. 

One can seek degree n >2 solutions within the ansatz (3.2) by replacing 9 with n9, as in 

[2], 

(p{r,9) = (-\/l - ^„(r)2cosn6', \/l - Zn{rY sin n9 , Zn{r)) (3.9) 

The profile function is then Zn{r) = — 1 + 2 exp(— r^/2n). But such fields are not even once 
differentiable at the origin, so are not genuine solutions of the Bogomol'nyi (or field) equation 
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in the sense that we demand. The problem is that tp has a conical singularity at (0, 0). To see 
this, let 

W^^^±^ (3.10) 

be the image of ip under stereographic projection from (0,0,-1). Note that 1^ is a good 
complex coordinate on a neighbourhood of (^(0, 0) — (0, 0, 1). Then, for this radially symmetric 
n-lump. 

Hence 

W^.(0,y) = yf(|jy \o{y') (3.12) 

which has a step discontinuity at y = 0. There is a similar problem with the radially symmetric 
degree n solutions obtained in [7]. Wc will sec below that genuine (at least twice differentiable) 
solutions of the Bogomol'nyi equation do exist for each n>l but constructing them requires 
some ingenuity. 

A geometric insight into the difficulty one faces can be obtained from Remark 2. In this 
case, the vacuum manifold is A^o — {(0,0,-1), so N' — N\Nq is a punctured sphere or, 
equivalently, an open disk. The deformed area form on N' is, in cylindrical coordinates, 

a = ^^ (3.13) 

1 + Z ^ ^ 

which gives A^' infinite total area. In fact, (A^', Q) can be visualized as a "cigar shaped" surface 
of revolution, with a single infinite cylindrical end replacing the missing point Aq, see figure 
1(a). This comes from identifying Q with the area form on the punctured sphere A^' associated 
with the metric 

where W is the stereographic coordinate defined in (3.10). The degree 1 energy minimizer 
constructed above can now be seen as an area-preserving diffeomorphism from to (A^', fl). 
The difficulty in constructing higher degree solutions is that any map of degree exceeding 1 
must have critical points. Any such critical point must get mapped to A'o, the end at infinity, 
and it is hard to arrange this while maintaining the area-preserving property of ip away from 
its critical points. Certainly ip cannot have any isolated critical points (as a generic map 
between 2-manifolds does), since we have the following proposition. 

Proposition 8 Let ip : M."^ S'^ be a solution of the model with U ((/?) = 1 + p^ satisfy- 
ing boundary condition (2.1). Then every connected component of the critical set of ip is 
unbounded. 

Proof: By Proposition 6, F'^ = f/^ everywhere and so p maps the critical set into 
Aq = U~^{0), the vacuum manifold. Assume, towards a contradiction, that has a bounded 
connected component Mi. Then for e > sufficiently small the closed 1-manifold (p'^^{—l + €) 
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(a) (b) (c) (d) 

Figure 1: The deformed target spaces {N',Q) embedded as surfaces of revolution, in the cases 
(a) t/ = 1 + (b) U=il + ^3)°•^ (c) U=il + ip.r-' and (d) U = (1 + ifsf-'il - ^sf-'- 

has a connected component F = S*^ whose interior contains Mi. Let S be the interior of F with 
Ml removed, and consider the restriction of to S. By Remark 2 this is an area-preserving 
surjective map from S to N,, = {ip : — 1 < (^93 < — 1 + e} with respect to fl. But S, being a 
bounded subset of M^, has finite area while {N^, Q) has inifnite area, a contradiction. □ 

Nonetheless, this model does have solutions in every homotopy class. We construct them 
as follows. Let ^ : (0, 00) x M — > be the diffeomorphism 

: {x,y) (-7- {\ogx,xy). (3.15) 

This map is area-preserving (with respect to the Euclidean metric on both spaces). Let 
: — 7- 5^ denote the unit lump solution constructed above, equations (3.2), (3.4). Then 
as remarked after Proposition 7, ip o ^ satisfies the Bogomol'nyi equation on the half-space 
(0, 00) X M. Clearly, lima;_>.o+ ■ip{£/{x,y)) = (0,0, —1) for all y. Hence the map 

.:M^-S^ .i.y) - (3.16) 

is continuous and satisfies the Bogomol'nyi equation away from the line x = 0. We claim that 
this is a genuine degree 1 solution of the Bogomol'nyi equation, and hence, the field equation. 
This amounts to the claim that (p is twice continuously differentiable everywhere. 

Proposition 9 The mapping y9 : — )■ S*^ defined in equation (3.16) is . 

Proof: Clearly, Lp is smooth away from the line x = 0, and all its derivatives vanish identically 
for X < 0. So it suffices to show that 

\Vx\, \Vy\, \Vxx\, \Vxy\, \Vyy\ (3-17) 
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all vanish in the limit x 0+, for all y. For a; > we have that (p{x,y) = iIj{X,Y) where 
X — log a;, Y = xy. Straightforward estimates using the explicit formulae (3. 2), (3. 4) yield 
that there exist constants C, i?* > such that for all R — -\/X^~+y^ > R^, 

\^xI\^y\ < CRe-^'l' 
\i1^xx\Mxy\Myy\ < CR\-''"'\ (3.18) 

Similarly, there exists constant X* < such that for all < a; < e^* , 



-X 



\X^\,\Xy\,\Y^\,\Yy\ < |y|+e 

l^xxl) l^xj/l) I ^2/2/ 1) l^^il) I ^^2/ 1) I ^^2/ 1 — ^ ■ (3.19) 

Hence, by the chain rule, for all < a; < = min{e"^*, e^^*} and all y 



< CRe-'^'/Wy] + e-^) < C{\X\ + ye^)e-^'/^(y + e"^) ^ (3.20) 



as a; — >■ 0"^, since then X — >■ — oo. Hence limj._^o+ \Vx{x-i y)\ — ^ for all y. The same argument 
deals with ^py. 

Turning to the second derivatives, we see from the chain rule and estimates (3.18), (3.19) 
that for all < a; < a;* and all y 



< C {i?e-«'/^e-2^ + R^e-'^'/^y^ + 6"^^) } 

< C(|X|+|/ + X2 + |/2)(^2^e-2^)e-^'/^^0 (3.2i: 



as a; — )■ 0"*", since then X — oo. Hence lim3,^o+ \Vxx{^-i y)\ — ^ for all y. The same argument 
deals with (pxy-i <Pyy n 

It seems likely that the mapping (f defined in (3.16) is actually smooth everywhere, but 
we have not proved this. Let us henceforth denote this degree 1 map, which satisfies 
the Bogomol'nyi equation everywhere, (^_|_. Note that E{ip^) = An, the topological minimum 
value in its homotopy class. Since it takes exactly the vacuum value on the left half-plane, one 
could call this solution a semi-compacton. However, by exploiting the invariance of E under 
area-preserving diffeomorphisms further, we can construct degree 1 energy minimizers with 
more tightly localized support. 

Consider the map 

: (0,oo) X (-^, ^) ^ (0,oo) X R, =e/'(a;, y) = (xcos^ y, tany). (3.22) 

Clearly is an area-preserving diffeomorphism. For any e > 0, denote by (p^^ the a;-translate 
of ip+ by e, that is, 

ip%{x,y) =(p+{x-e,y). (3.23) 

The support of ip^ is the half plane x > e. Denote by 5" the infinite half-strip S = (0, oo) x 
(— |, I), and consider the mapping 
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Figure 2: The energy density of a semi-compacton. 



By construction, this is continuous everywhere and on the complement of dS, the boundary 
of the strip S. It also satisfies the Bogomol'nyi equation on '^^\dS. By construction, its 
support is a subset of the closure of S*. In fact. 



Hence is constant on a neighbourhood of OS, and so is trivially on OS. Hence is 
everywhere. By construction, E{ip^) = 4tt, that is, is a degree 1 energy minimizer, 
which we call a semi-compacton. It has a single energy density maximum located at the point 
(1 + e, 0). The energy density along the line y = is 



So has an energy tail which decays along the strip 5* like x faster than any power, 
but slower than exponential. The energy density of (for e very small) is plotted in figure 
2. 

By precomposing (f^ with an area preserving diffeomorphism 



where a > 0, we can construct semi-compactons with support in an arbitrarily thin, half 
infinite strip. Similarly, the strip can be deformed to follow any non-self-intersecting half 
infinite curve which escapes to infinity. The mapping '^□(x,?/) = '/'□(a;, — ?/) is an anti-semi- 
compacton, of degree —1. By gluing together (anti-)semi-compactons with disjoint support, 
one obtains energy minimizers in every homotopy class. Gluing together > semi- 
compactons and n_ > anti-semicompactons yields degree n = n+ — n_ fields which, by 
Proposition 5, are solutions of the field equation, but have energy 4(n+ -|- n_)TT > 4|n|7r. 
So each homotopy class contains critical points of E of arbitrarily high energy. Even more 
surprising, these critical points are not saddle points of E but are, in a certain sense, linearly 
stable. 

To see this, one must construct the Hessian operator for the functional E{ip) based at a 
critical point if. We recall that this is defined as follows. Let cps^t be a two-parameter variation 
of a critical point ip : M ^ N of E, and let X = dsips,t\s=t=o, Y = dt^s,t\s=t=o G T{lp~^TN) be 



supp(y9Q = {(x, ?/) : X > e/ cos^y} C 5". 



(3.25) 




(3.26) 



£^ (x, y) = {ax, a y) 



(3.27) 
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the associated inifnitesimal variations. Then the Hessian oi E ai ip is the symmetric bihnear 
form 

Hess(X,r) = — 



dsdt 



(3.28) 

s=t=0 



on r((/9 ^TN). The associated Hessian operator is the self-adjoint hnear differential operator 
Jif : T{ip-^TN) T{ip-^TN) defined such that 



Hess(X,F)= / h{X,jrY)vo\M = {X,jrY). (3.29) 
Jm 

One uses the spectrum of Jif to classify the critical point cp. In particular, if Jif has both 
negative and positive eigenvalues, is a saddle point. If the quadratic form Hess(X, X) is 
non-negative, one says that ip is linearly stable (although (p may actually be dynamically 
unstable; e.g. is a linearly stable critical point of /(,x) = —x'^). 
For the energy under consideration here, one finds that [10] 

jfy = - J(V|^r + dip{fi5dip*iY0j)) + (V^ grad V) o ip (3.30) 

where = ^dip*LJ G T{TM), is the Levi-Civita connexion on TN, is its puUback to 
(p-'^TN, V — |C/^ and l denotes interior product {laco — co{A, •)). The exact details of this 
formula are not important. We will need only the following Lemma. 

Lemma 10 Let ip : M ^ N be a critical point of E, Jif be its Hessian operator and x G 
(M\supp(^). Then, for all Y e r{(p-'^TN), 

{J^Y){x) = 0. 

That is, the Hessian operator vanishes identically off the support of (p. 

Proof: The complement of supp ip is open by definition, so ip is constant on a neighbourhood 
of x. It follows that — and d(/7 = on a neighbourhood of x, so the first two terms in 
J^fY vanish at x for all Y. Consider now the zeroth order piece 

J^y = (Vy grad V) o ip. (3.31) 

Since V ^ ^U^, 

Jif^Y = {Vy{U grad U)) o ip = {h{Y, grad U) grad U) o ip + {UVy grad U) o ip. (3.32) 

But U{ip{x)) = (since x ^ supp(/9) and, since U is assumed non-negative, ip{x) is a minimum 
of U, and hence (grad t/) (99(2;)) = also. Hence, for all Y, {MqY){x) = □ 



Proposition 11 Let ip : M ^ N be any solution of (2.14) constructed by superposing (anti-) 
semi-compactons ipi with support in disjoint strips Si, i = 1,2, ... ,m. Then for all Y G 

r-\TN), 

Hess(y, Y) > 0. 
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Proof: By Lemma 10, 

m „ 

Hess(r,r) = V/ h{Y,J^Y)vo[M (3.33) 

i=l J Si 

where ^ denotes the Hessian operator associated to the (anti-)semi-compacton Each 
term in this sum is non-negative for all Y . For if not, then, by Lemma 10, there exists i and 
a section Y such that 

/ j^r)voiM = / J^r)voiM < o, (3.34) 

J M J Si 

which contradicts the fact that ipi minimizes E in its homotopy class. □ 



Physically, the point is that semi-compactons exert no forces on one another, so 'Piz'fic: 
superpositions are (marginally) stable, by stability of their constituent parts. 

So this model supports degree n (marginally) stable multi-semi-compactons of energy 
4(|n| + 2k)7i for all n G Z and k G Z>o. All these solutions have (multiple) tails escaping 
to infinity, along which the energy density decays like x~^°^^. 



4 Compactons revisited 

When will an extreme baby Skyrme model support genuine compactons? The geometric 
picture outlined above immediately gives a necessary condition on U, namely that N' = N\No 
(the target space with its vacua removed) should have finite volume with respect to the 
deformed area form fl = uj/U. 

Proposition 12 Let : M.^ ^ N be a surjective solution of T{(f) — of compact support. 
Then {N', fl) has finite volume. 

Proof: By Proposition 6, — (U o (p^ everywhere. Since F^{x) = if and only if a; G M^, 
(f defines an area preserving map from each connected component of M\M^ into (N', Q). The 
union of the ranges of each such map is all N' (since (p is surjective), and hence the area of 
N' cannot exceed the area of M\M^ C supp ip. □ 

Conversely, if {N', fi) has finite area A, let M' be any subset of M = of area A which 
is diffeomorphic to N'. For example, if Nq consists of p vacua, one could take M' to be an 
open disk of area A + pe with p small disjoint closed disks of area e removed. Construct an 
area-preserving diffeomorphism tjj : M'\N', using the method of Moser, for example [6], and 
extend t/j to the whole of M by a piecewise constant map on M\M'. This map cp certainly 
has compact support, and satisfies the field equation except, perhaps, on the boundary of M'. 
Hence (/? is a genuine solution if and only if it is C^. 

For example, consider the model with 

U{<p)^{l + cpsr (4.1) 
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where ^ < a < 1. Here A^' = 5'^\{(0,0, —1)}, diffeomorphic to an open disk, but, unlike the 
case a — 1 considered in the previous section, N' now has finite volume. 

It can be visualized as a baloon shaped surface of revolution, with a conical singularity at the 
missing vacuum point, see figure 1(b), (c). An obvious choice for the open set M' is the disk 
of radius R = 2i-°/2(l - a)- 2. There is an area-preserving diffeomorphism M' — >■ N' within 
the radial ansatz (3.2), 



z{r) 



2'-'^- 1(1 -ay 



1 

l-a 



- 1, (4.3) 



which, when extended by (0,0,-1) outside the disk M' gives a map — )■ 5'^ of degree 
1 solving the field equation everywhere. This (up to reparametrization) is the compacton 
reported by Adam et al [2]. Note that one can obtain its key qualitative features without 
solving any equations, e.g. it occupies area jz^tt and has total energy 

Another interesting choice is 

U{^) = {l + ^^r{l-^,f (4.5) 

where a,/? e [|, !)• Now N' is diffeomorphic to a cylinder and has finite total area A(q;,/3), 
a complicated function of a, j3 involving hyper geometric functions. An embedding of A^' as a 
surface of revolution in the case a — 0.5, ^ — 0.7 is depicted in figure 1(d). One can take M' 
to be any annulus of total area A, 

M' = {(x, y) : Rl<x'^ + y'^ < Rl} where n{Rl - Rj) = A, (4.6) 

and construct an area-preserving diffeomorphism ip : M' — )■ A^' within the ansatz (3.2), then 
extend this by (0,0, —1) for r < Ri, and (0,0, 1) for r > i?2- It is straightforward to check 
that this field is C^, and hence defines a ringlike compacton. By choosing R2 sufficiently large 
(and Rl close to R2) this ring can be arbitrarily big. Hence one can construct n-compactons, 
with n rings nested inside one another, as well as the more obvious multi-ring solutions. Adam 
et al consider only the degenerate case that the annulus is a punctured disk [2] , so we shall go 
through this construction in more detail. 

The deformed area form (in cyhndrical coordinates) is Q = (1 -|- Z)-"(l — Z)~^dQ A dZ, 
so a field within the ansatz (3.2) satisfies the Bogomol'nyi equation if and only if 

-r. (4.7) 



{l + zY{l-zy 
Define the function 



Q : [-1, 1] ^ [0,A/{2n)], Q{Z) = £ ^ _ . (4.8) 
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Then Q ^ is an increasing, surjective map [0,/l/(27r)] [—1, 1], and 



z{r) = Q-\C -''^) (4.9) 

solves (4.7) for any constant C > 0. We require (/9(0,0) = (0,0, 1), so insist that C > A/(27r). 
Set Ri = \pKi and define i?2 > -^i such that 7r(i?| — i?f) = A. Then the extended profile 
function is 

{1 < r < i?i 

Q-W{R\-r'')) Ri<r<R2 (4.10) 
-1 r > i?2 

Note that the associated field has support in an annulus of total area A, as expected. Note 
also that it is constant in a neighbourhood of the (polar) coordinate singularity at r = 0, so 
to check that ip is C^, it suffices to check that z{r) is C^. This is clear, except at the points 
r — Ri, where z — —1 and r = R2, where z = 1. By the Bogomol'nyi equation, 

z'{r) = -r(l + z(r))"(l - z(r))^ (4.11) 

on (i?i,i?2), whence \im^^j^+ z' (r) = lim^^j^- z'{r) = 0. Hence z{r) is C^. Differentiating 
(4.11), 

z"{r) = -{1 + z{r))''{l - zir))" 

-r^{a{l + z{r)f''-\l - z{r)f^ - /3{1 + z{r)f''{l - z{r)f^-^} (4.12) 

on (i?i,i?2), whence lim^_^^+ z"{r) — lim^_^^- z"{r) — also (note a,/3 > |). Hence z{r) is 

One can precompose this map with an arbitrary area-preserving diffeomorphism : — )■ 
to obtained deformed ring-like compactons. Choosing R2 very close to Ri, then deforming, 
produces closed string-like compactons. In fact, one can precompose it with a degree n area- 
preserving covering map 

^ : M2\|(o,0)} ^ i?2\|(0,0)}, ^ : {r,e) ^ {n-'^r,ne) (4.13) 

to obtain a degree n annular compacton. Unlike the single vacuum case, this is still (even 
at the origin) because the degree 1 compacton is constant on a neighbourhood of the origin. 

Finally, consider the case of potential (4.5) in the case /3 = 1. This supports a de- 
gree 1 energy minimizer which decays like exp(— r^/2) to (0,0, 1) as r ^ 00, and is exactly 
(0,0, —1) on any closed disk centred on the origin. By precomposing this with appropriate 
area preserving maps, as in the previous section, we can produce a semi-compacton localized 
in a semi-infinite strip, with a x~^°^^ tail, but with a hole (of any finite area) in the middle 
of the lump, where it has exactly zero energy. Clearly, by introducing more vacua, one can 
dream up models with even more bizarre energy minmizers. 
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5 Concluding remarks 



We have shown that the extreme baby-Skyrme model with energy 

E^l f {b • X ify)]^ + U{ipf}dxdy, U{ip) = l + ips (5.1) 

supports, in every homotopy class, semi-compacton solutions of quantized energy E = 47r(|n| + 
2k) where n e Z is the degree of (/? and A; is a non-negative integer. These solutions are (at 
least) twice continuously differentiable everywhere, and consist of |n| + 2k (anti-)lumps, each 
localized in a semi-infinite strip. Each lump has a tail escaping to infinity, along which the 
energy density decays like s^^°s*, where s is a length variable along the strip. All these 
solutions are at least marginally stable, and when A; = 0, are global energy minimizers in their 
homotopy class. 

Replacing the potential term hy U — (l + f/^a)", | < a < Iwe have given a geometric 
interpretation to the construction of compactons proposed in [2] , and clarified the conditions 
under which these are (hence classical solutions of the field equation). In the case of two- 
vacuum potentials U = {1 + (^3)^(1 — ^3^, we have constructed annular compactons, and 
described how these can be embedded inside one another, and deformed into closed string-hke 
solutions. 

It is interesting to compare this situation with the case where M is compact. The role 
of the potential term ^f/^ on is to prevent lumps dissipating by spreading indefinitely. 
On compact M, the very compactness of M docs this job, so one might expect that similar 
results (existence of minimizers in every homotopy class) might hold here in the simple case 
U — 0. This turns out to be entirely false. Indeed, it was shown in [11] that all critical 
points of Ei[ip) on a compact Riemann surface have ip*u: coclosed. Now ip*u is automatically 
closed for all 93 (since (i(p*uj = (p*duj = 0), so if 9? solves the field equation for E4, ip*cu is 
harmonic. Hence, by the Hodge Theorem, ip*u} = constant x voIm, that is, (f : M ^ N 
is, up to a homothety of {M,g), an area-preserving covering map, or £'4(93) = 0. So if the 
target is N — S^, any solution either has degree 0, or is an area-preserving diffeomorphism 
M ^ S"^ (since S'^ is simply connected, any covering map is a diffeomorphism). It follows that 
if M = S"^, the model has solutions only in the degree —1, 0, 1 classes, while if M is any other 
compact Riemann surface, it has only trivial (degree 0, energy 0) solutions. The contrast with 
M = and U is striking. 

The results of this paper raise two obvious interesting questions. First, can one understand 
the moduli space of degree 1 energy minimizers of this model? What about the reduced 
moduli space, that is, the set of minimizers modulo the action of the group of area-preserving 
diffeomorphisms of M^? Clearly, the radially symmetric lump ip, the half lump and the 
semi-compacton are three different points in this space. Do they lie in the same connected 
component? Is the moduli space, in fact, connected? If so, can it be given a manifold structure? 
If not, can its components be enumerated? Such questions are mathematically well-defined 
(for example, we can give the set of all maps the compact-open topology, the moduli space the 
relative topology from this, and the reduced moduli space the quotient topology from this) 
but seem formidably challenging. 

Second, can one study the dynamics of semi-compactons? This question is rather subtle. 
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because the Euler-Lagrange equation descending from the obvious Lorentz-invariant time- 
dependent extension of the model, with Lagrangian density 

^^\Vp- {d,^ X dM\[^ ■ {d^^ X 9»] - \u{^f (5.2) 

is not a true cvohition equation. The problem is that, at any spatial point {x,y) G M? where 
ipx-i do not span T^S^ (that is, at any critical point of (p(t,-) : S^), the fields ip 

and (fit, do not uniquely determine (pu- In particular, the Cauchy problem for any initial data 
</?(0), (pt{0) is ill-defined if (p{0) has any critical points. This is immediately a problem for any 
initial data of degree > 2, since any such field has critical points by topological considerations. 
For semi-compactons, the problem is particularly severe, since these are critical on unbounded 
regions of M^. If the moduli space of semi-compactons can be understood, one could perhaps 
study the dynamics of a single semi-compacton within the geodesic approximation. There 
are some indications that the kinetic energy functional of (5.2) equips the moduli space, at 
least formally, with an incomplete Riemannian metric. Less speculatively, one could abandon 
Lorentz invariance (which is, in any case, an unnatural assumption for condensed matter 
applications) and give the model the usual kinetic energy term, that is, 

^=l^f^t- ■ X ^,)]2 - 1^/(^)2. (5.3) 

The Euler-Langrangc equation is now a genuine evolution equation, although it is not tech- 
nically hyperbolic. It would be interesting, and numerically straightforward, to study the 
scattering of semi-compactons in this model. 
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